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Prediction of inception of sheet cavitation on solid walls has been recognized to be very difficult, since it is sig- 
nificantly affected by the boundary layer flow characteristics, the population of free nuclei, the nuclei held in the 
wall roughness, the amount of dissolved air in liquid and so on. It has not sufficiently been made clear how the 
inception is affected by the conditions of water qualities and background flow characteristics. In this study, high 
speed observation of inception of sheet cavity from free nuclei is conducted for a two-dimensional conver- 
gent-divergent nozzle flow, where the sheet cavity forms just downstream of the nozzle throat. The effects of the 
amount of dissolved air and the free stream velocity on the inception process of sheet cavitation is examined. In 
addition, the bubble nuclei density, which is well known to be important factor for cavitation inception, is pas- 
sively controlled by the filter installed in the tunnel. From the observations, it is confirmed that the nuclei number 
density significantly affects the formation of sheet cavity rather than the other two parameters. In conditions with 
large nuclei number density, the sheet cavity does not form, and bubbly cavitation appears instead. In the case 
with small nuclei number density, the sheet cavity forms from a single flowing nucleus and develops streamwise- 
ly and spanwisely. In the conditions with medium nuclei number density, the sheet cavity also forms but is short- 


er/narrower streamwisely/spanwisely, due to interaction of other nuclei flowing near the formed sheet cavity. 
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layer flow if the cavitation near solid surfaces is focused 
on. The cavitation inception near solid surfaces has been 
studied by many researchers (for example [2-4]). How- 
ever, about the inception of sheet cavitation, there is still 
a review article published [5], and it seems to be a com- 


Introduction 


Cavitation is a phase change phenomenon from liquid 
to gas phases due to the hydrodynamically reduced pres- 
sure below its saturated vapor pressure level [1]. It is 


well known that, in addition to the low pressure condition, 
the inception of cavitation requires the bubble nuclei 
from which the vapor bubbles explode as well as the suf- 
ficient time which allows the nuclei to grow. The former 
one, the nuclei distribution, is dependent on the flow 
system such as liquid (water) quality including the 
amount of dissolved air, trapped air in the roughness 
element of solid surface, and so on. The latter is well 
affected by the flow characteristics, especially boundary 
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mon understanding that there exist many scenario for the 
inception of sheet cavitation but none of which is the 
typical one, which probably depends upon the flow cha- 
racteristics and the water quality as we stated before. 

In this study, high speed observation of inception of 
sheet cavity from free nuclei is conducted for a two- 
dimensional convergent-divergent nozzle flow, where the 
sheet cavity forms just downstream of the nozzle throat. 
The effects of the amount of dissolved air and the free 
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Nomenclature 
h height of nozzle throat 
p pressure (Pa) 
Py saturation pressure (Pa) 
U velocity (m/s) 
Greek letters 
p density (kg/m°) 


stream velocity on the inception process of sheet cavita- 
tion is examined. In addition, the bubble nuclei density, 
which is well known to be the important factor for cavi- 
tation inception, is passively controlled by the filter in- 
stalled in the tunnel. By using the image analysis of 
flowing nuclei, the nuclei density distribution is obtained. 
Then the inception process of sheet cavitation is dis- 
cussed. 


Experimental Procedure 


Experimental apparatus 


Figure 1 shows our small cavitation tunnel used for 
this study, which consists of an upstream tank, a test sec- 
tion, two circulating pumps operated in parallel and a 
turbine flow meter. There are two lines split after the tur- 
bine flow meter; the flow in the one line directly returns 
to the tank via the flow control valve, while the flow in 
the other does via the air removal tank with wounded 
polypropylene cartridge filter, whose nominal pore di- 
ameter is lum. The upstream tank is connected to the 
atmosphere via a vacuum pump at the top so that the 
system pressure can be adjusted to the desired pressure. 

Figure 2 shows schematic view of the test section. The 
test section has a rectangular cross-section with a height 
of 20 mm and a width of 20 mm, and the convergent- 
divergent nozzle is installed which consists of a top 
straight wall and a bottom inclined wall with the diver- 
gent angle of 8.4 degrees. The nozzle throat height is 10 
mm, and the length of throat portion is 3 mm. Sheet cav- 
ity initiates slightly downstream of the nozzle throat and 
develops along the bottom nozzle wall. All walls are 
made of acrylic resin, which enables us to set cameras 
and lightings from all directions. A pressure transducer is 
installed at UP of the top wall to measure the working 
pressure pyp of the flow. 


Experimental procedure 


In the experiment, the working pressure is gradually 
decreased from non-cavitation condition, and the cavita- 
tion inception from flowing nuclei is observed from top 
and side by two high-speed cameras. 

For the top view, the nozzle throat region indicated by 
area enclosed by dotted lines in Fig.2 is recorded with a 
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o cavitation number 
Subscripts 

l liquid phase 

th nozzle throat 

UP nozzle inlet 


back lighting method by Vision Research, Phantom V310 
with an exposure time of lus. The spatial resolution is set 
to be 1280 x 400 pixels, resulting in the actual resolution 
of 16 um/pixel. From the recorded movie, the diameter of 
each bubble nucleus is measured by an image analysis 
described later. On the other hand, for the side view, the 
region indicated by area enclosed by dashed lines in Fig. 
2 is recorded with a back lighting method by Vision Re- 
search, Phantom V4.3 with an exposure time of 10 us. 
The spatial resolution is 608 x 152 pixels, which corres- 
ponds to the actual resolution of 12 um/pixel. These two 
observations are made simultaneously with the frame rate 
of 6259 fps, which enables us to trace the motion and 
growth of each nucleus. 

In this experiment, the area averaged velocity at noz- 
zle throat is set to be Up= 10 or 4m/s. Corresponding 
Reynolds number based on the nozzle throat velocity Un 


Air supply 


Vacuum pump 


Test section 


Circulating pumps 


Fig. 1 Experimental apparatus. 
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Fig. 2 Test section and recorded area with high speed cam- 
eras. 
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and the nozzle height h is roughly 10° and 2 x 10* re- 
spectively. Before the experiment, the amount of dis- 
solved air in water is roughly controlled by deaeration, 
and Dissolved Oxygen (DO) in water is measured before 
and after the experiment. In some conditions of DO and 
Un, the downstream line is switched to those with or 
without air removal filter, in order to see the effects of 
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bubble nuclei distribution on the inception of sheet cavity. 


By using the pressure at UP, pup, and Um, the cavitation 
number is defined as 


g2 Pup — Py 
z 2 
PU n / 2 
where p, is a saturated vapor pressure. 


(1) 


Image processing 


To measure the diameter of bubble nuclei, the image 
analysis is carried out in the way shown in Fig.3. The 
images from the recorded movie in top view () are used. 
To observe only the nozzle throat region, the other area is 
masked as shown in @. All recorded images, about 
3,000 frames with the time duration of 0.5s, are averaged 
to construct the background image ©. Then, the back- 
ground image is subtracted from all original image to 
obtain the images @) which contain only the information 
of bubble nuclei. Edge detection is carried out with some 
threshold value © and the enclosed area is filled with 
white nuclei ©. Finally, the bubble nuclei located on the 
boundaries are removed ©. 

From the above-obtained images, we calculate the 
projected area of bubble nuclei by counting the pixels 
occupied by the bubble image. Then, the radius of each 
bubble nucleus is converted from the projected area un- 
der the assumption of spherical bubble nuclei. When the 
sufficient amount of the bubble nuclei data is obtained, 
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we can calculate the nuclei density distribution against 
the assumed reference volume of 3 mm (length) x 20 mm 
(width) x 10 mm (depth). It should be noted that the 
measured radius of nuclei is larger in some extent than 
the upstream one due to the flow acceleration from up- 
stream to the measured location. 


Results and Discussion 


Characteristics of nuclei density distributions 


Figure 4 shows the number density distribution func- 
tion of bubble nuclei in cases without using the air re- 
moval filter. The cavitation number is set to be around 
the inception point of sheet cavity. Three different DO 
conditions, high (DO = 8-9 mg/L, almost saturated at the 
atmospheric pressure), medium (DO = 4-5 mg/L, par- 
tially degassed) and low (DO = 1-2 mg/L, well degassed) 
conditions, are examined. The nozzle throat velocity is 
basically set to be Un= 10 m/s but for low DO condition, 
the velocity condition of Un= 4 m/s is also examined. In 
terms of dissolution of gas in liquid, the relation between 
the dimensional working pressure and the saturated 
pressure determined by Henry’s law is important. Only in 
the low DO condition with high nozzle throat velocity 
(plotted by reverse triangle in the figure), the dissolved 
air is under-saturated, while in the other three conditions 
over-saturated. Actually, as can be seen in the figure, the 
number of bubble nuclei is so small in the low DO condi- 
tion with high nozzle throat velocity. The number of 
bubble nuclei is apparently insufficient to obtain the re- 
liable nuclei number density function; longer recording 
duration is necessary to do so. On the other hand, the 
number density of bubble nuclei is larger over all bubble 
radius regime for the other three cases and it is the largest 
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Fig. 3 Procedure of image processing for measurement of radius of bubble nuclei 
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for the high DO condition (DO = 8.4 mg/L). In this con- 
dition, the over-saturation (working pressure minus Hen- 
ry’s pressure) is the largest, which results in the largest 
number density of bubble nuclei. The exponent of the 
number density distribution to the radius of nuclei is 
roughly -2.1 to -2.8, which is consistent to the published 
data (for example, [1]). 

Figure 5 shows the number density distribution func- 
tion of bubble nuclei in cases with the air removal filter 
used. The number density distribution of bubble nuclei is 
similar to the previous one in the under-saturated condi- 
tion, that is the low DO condition with higher nozzle 
throat velocity U;,= 10 m/s. On the other hand, it is no- 
ticed that the number density distributions for three 
over-saturated conditions coincidently agree well with 
each other. The exponent of the number density distribu- 
tion to the radius of nuclei is roughly to -1.8. It is be- 
lieved that, by the air removal filter, large bubbles can be 
removed but smaller ones whose diameter is less than the 
order of 1 um (pore diameter of the filter) can pass 
through the filter (or possibly generated by the filter it- 
self). In over-saturated conditions, those tiny bubbles 
might be able to grow by the mass diffusion of the dis- 
solved air through the surface of nuclei, which could be 
the reason for the coincidence of the nuclei density dis- 
tribution. 

Figure 6 shows the number density of bubble nuclei, 
that is the number of nuclei in unit volume, plotted 
against cavitation number o. It is seen that the data can be 
separated roughly into three groups, depending upon the 
nuclei density distribution; the large nuclei density con- 
dition which corresponds only to unfiltered condition 
with high DO (open square), the small nuclei density 
conditions which correspond to the filtered/unfiltered 
conditions of low DO with high nozzle throat velocity 
Umn= 10 m/s (closed and open reverse triangles) and the 
medium nuclei density conditions which correspond to 
the other conditions, in which the nuclei density distribu- 
tion functions coincide with each other as has been seen 
in Figs. 4 and 5. Since the plotted range of cavitation 
number is the range of the inception of sheet cavity ex- 
cept for filtered condition with high DO as described 
later, the nuclei density seems to be the primary factor for 
the inception of sheet cavitation. 


Inception behavior of sheet cavity 


Figure 7 shows the typical images of bubble nuclei 
behavior for (a) non-filtered and (b) filtered conditions 
with high DO (= 8.4 mg/L). In the non-filtered condition, 
there can be seen so many bubble nuclei. The sheet cavi- 
ty is not at all formed, and the bubbly cavitation appears 
instead. On the other hand, in the filtered condition 
shown in Fig. 7 (b), there can be seen apparently a few 
nuclei. A nucleus travels the distance of 1.6 mm between 
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0 ms to 0.16 ms along the wall, and its size is approx- 
imately doubled. Then the nucleus attaches to the wall 
downstream of the throat, and forms finger-like sheet 
cavity at 9.7 ms. It splits and grows in streamwise and 
spanwise directions. 
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Fig. 4 Density distribution of bubble nuclei number for non- 
filtered cases. 
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Fig. 5 Density distribution of bubble nuclei number for fil- 
tered cases. 
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Fig.6 Bubble nuclei number per unit volume. 
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Figure 8 shows the inception process of sheet cavita- 
tion for the filtered middle DO condition (DO = 4.3 mg/L) 
with the nozzle throat velocity of Up= 10 m/s. The in- 
ception process in this middle DO case is basically simi- 
lar to that for the filtered high DO condition shown in Fig. 
7(a). A few flowing nuclei is observed in those two con- 
ditions, and the sheet cavity starts to form from a single 
nucleus. It is observed that the sheet cavity grows but its 
extent is limited in both the streamwise and spanwise 
directions. When the other nuclei pass over the sheet cav- 
ity, a part of the sheet cavity is torn off as shown in Fig. 8 
(b). Although the DO conditions are different, the 
processes of sheet cavity inception are similar when the 
number density distributions of nuclei are similar, and the 
growth of the sheet cavity is limited by the interaction 
of .flowing nuclei. 

Figure 9 shows the typical pictures of inception of 
sheet cavitation for the non-filtered low DO conditions 
with the nozzle throat velocity of (a) U,,= 10 m/s and (b) 
Ui, = 4 m/s. In both cases, the sheet cavity is seen to form 
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(a) Non-filtered condition (ø = 0.85) 
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from a single flowing nucleus. In the case with Up= 10 
m/s, a nucleus travels the distance of 1.8 mm along the 
wall in a time duration of 0.16 ms. It grows faster com- 
pared with the high DO condition as has been shown in 
Fig. 7 (b). After the nucleus attaches to the wall, finger- 
like sheet cavity grows more widely in streamwise and 
spanwise directions than the high middle DO conditions, 
probably due to less number of flowing nuclei. In the 
case with U;,= 4 m/s, a spherical nucleus travels the dis- 
tance of 2.8 mm along the wall in a time duration of 0.64 
ms. It grows more rapidly than that observed for U;,= 10 
m/s, which seems to be simply because the residence 
time of the bubble nucleus in low pressure region is 
longer for the lower velocity case. The tail of the nucleus 
attaches to the wall, and the sheet cavity forms. 

It has been observed that the moving velocity of bub- 
ble nuclei is a little bit faster than the nozzle throat veloc- 
ity Um which seems to indicate that there exists a slip 
velocity between the bubble nuclei and the surround liq- 
uid, even though the bubble nuclei are small. For more 
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Fig. 7 Comparison of bubble nuclei behaviors in high DO condition (Um = 10 m/s, DO = 8.4 mg/L) 
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(a) Inception from a free nucleus 
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Fig.8 Sheet cavitation inception in filtered middle DO condition (Um = 10 m/s, o = 0.86, DO = 4.3 mg/L) 
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(a) Non-filtered condition with Un = 10 m/s 
(o = 0.82, DO = 1.6 mg/L) 
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(b) Non-filtered condition with Un = 4 m/s 
(o = 0.82, DO = 1.4 mg/L) 


Fig.9 Effect of nozzle throat velocity on inception of sheet cavitation (Low DO condition) 


meaningful discussion about the slip velocity, it is neces- 
sary to specify the local velocity of liquid phase, which 
remains in our future study. 


Conclusions 


In this study, visual observations of inception of sheet 
cavitation are carried out for a two-dimensional conver- 
gent-divergent nozzle flow. The effects of the amount of 
dissolved air, the nuclei number density, and the flow 
velocity is investigated. It is found that the nuclei number 
density is the primary factor for the inception of sheet 
cavity from free nuclei. 

In large nuclei number density condition, the sheet 
cavity does not form and the bubbly cavitation appears 
instead. Many bubble nuclei interact with each other, 
which prevents the formation of sheet cavity from each 
nucleus. On the other hand in the low nuclei number 
density condition, the sheet cavity starts to form from a 
single nucleus and it grows both streamwisely and span- 
wisely. In the middle nuclei number density condition, 
the sheet cavity starts to form from a single nucleus re- 
gardless of the amount of dissolved air. Its extent in 
streamwise and spanwise directions is limited due to the 
interaction with other flowing nucleus, perhaps depend- 
ing upon the local/instantaneous bubble nuclei distribution. 
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